Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids.

Discovery and analysis of genetic variants underlying agriculturally important traits are key to molecular breeding of crops. Reduced representation approaches have provided cost-efficient genotyping using next-generation sequencing. However, accurate genotype calling from next-generation sequencing data is challenging, particularly in polyploid species due to their genome complexity. Recently developed Bayesian statistical methods implemented in available software packages, polyRAD, EBG, and updog, incorporate error rates and population parameters to accurately estimate allelic dosage across any ploidy. We used empirical and simulated data to evaluate the three Bayesian algorithms and demonstrated their impact on the power of genome-wide association study (GWAS) analysis and the accuracy of genomic prediction. We further incorporated uncertainty in allelic dosage estimation by testing continuous genotype calls and comparing their performance to discrete genotypes in GWAS and genomic prediction. We tested the genotype-calling methods using data from two autotetraploid species, Miscanthus sacchariflorus and Vaccinium corymbosum, and performed GWAS and genomic prediction. In the empirical study, the tested Bayesian genotype-calling algorithms differed in their downstream effects on GWAS and genomic prediction, with some showing advantages over others. Through subsequent simulation studies, we observed that at low read depth, polyRAD was advantageous in its effect on GWAS power and limit of false positives. Additionally, we found that continuous genotypes increased the accuracy of genomic prediction, by reducing genotyping error, particularly at low sequencing depth. Our results indicate that by using the Bayesian algorithm implemented in polyRAD and continuous genotypes, we can accurately and cost-efficiently implement GWAS and genomic prediction in polyploid crops.

Incidence of RNA viruses infecting taro and tannia in East Africa and molecular characterisation of dasheen mosaic virus isolates.

Taro (Colocasia esculenta) and tannia (Xanthosoma sp.) plants growing in 25 districts across Ethiopia, Kenya, Tanzania and Uganda were surveyed for four RNA viruses. Leaf samples from 392 plants were tested for cucumber mosaic virus (CMV), dasheen mosaic virus (DsMV), taro vein chlorosis virus (TaVCV) and Colocasia bobone disease-associated virus (CBDaV) by RT-PCR. No samples tested positive for TaVCV or CBDaV, while CMV was only detected in three tannia samples with mosaic symptoms from Uganda. DsMV was detected in 40 samples, including 36 out of 171 from Ethiopia, one out of 94 from Uganda and three out of 41 from Tanzania, while none of the 86 samples from Kenya tested positive for any of the four viruses. The complete genomes of nine DsMV isolates from East Africa were cloned and sequenced. Phylogenetic analyses based on the amino acid sequence of the DsMV CP-coding region revealed two distinct clades. Isolates from Ethiopia were distributed in both clades, while samples from Uganda and Tanzania belong to different clades. Seven possible recombination events were identified from the analysis carried out on the available 15 full-length DsMV isolates. Nucleotide substitution ratio analysis revealed that all the DsMV genes are under strong negative selection pressure.

Molecular characterisation of a putative new polerovirus infecting pumpkin (Cucurbita pepo) in Kenya.

Next-generation sequencing of RNA extracted from a pumpkin plant with mosaic symptoms in Kenya identified the presence of a polerovirus sequence closely related to pepo aphid-borne yellows virus (PABYV). The near-complete polerovirus sequence comprised 5,810 nucleotides and contained seven putative open reading frames (ORFs) with a genome organisation typical of poleroviruses. BLASTp analysis of the translated sequences of ORFs 0, 1 and 2 revealed that their amino acid sequences differed by more than 10% from the corresponding protein sequences of other poleroviruses. These results suggest that this virus is a putative novel member of the genus Polerovirus, which has been provisionally named "pumpkin polerovirus" (PuPV).

Viral metagenomics of aphids present in bean and maize plots on mixed-use farms in Kenya reveals the presence of three dicistroviruses including a novel Big Sioux River virus-like dicistrovirus.

BACKGROUND: Aphids are major vectors of plant viruses. Common bean (Phaseolus vulgaris L.) and maize (Zea mays L.) are important crops that are vulnerable to aphid herbivory and aphid-transmitted viruses. In East and Central Africa, common bean is frequently intercropped by smallholder farmers to provide fixed nitrogen for cultivation of starch crops such as maize. We used a PCR-based technique to identify aphids prevalent in smallholder bean farms and next generation sequencing shotgun metagenomics to examine the diversity of viruses present in aphids and in maize leaf samples. Samples were collected from farms in Kenya in a range of agro-ecological zones. RESULTS: Cytochrome oxidase 1 (CO1) gene sequencing showed that Aphis fabae was the sole aphid species present in bean plots in the farms visited. Sequencing of total RNA from aphids using the Illumina platform detected three dicistroviruses. Maize leaf RNA was also analysed. Identification of Aphid lethal paralysis virus (ALPV), Rhopalosiphum padi virus (RhPV), and a novel Big Sioux River virus (BSRV)-like dicistrovirus in aphid and maize samples was confirmed using reverse transcription-polymerase chain reactions and sequencing of amplified DNA products. Phylogenetic, nucleotide and protein sequence analyses of eight ALPV genomes revealed evidence of intra-species recombination, with the data suggesting there may be two ALPV lineages. Analysis of BSRV-like virus genomic RNA sequences revealed features that are consistent with other dicistroviruses and that it is phylogenetically closely related to dicistroviruses of the genus Cripavirus. CONCLUSIONS: The discovery of ALPV and RhPV in aphids and maize further demonstrates the broad occurrence of these dicistroviruses. Dicistroviruses are remarkable in that they use plants as reservoirs that facilitate infection of their insect replicative hosts, such as aphids. This is the first report of these viruses being isolated from either organism. The BSRV-like sequences represent a potentially novel dicistrovirus infecting A. fabae.

Analyses of mitochondrial genes reveal two sympatric but genetically divergent lineages of Rhipicephalus appendiculatus in Kenya.

BACKGROUND: The ixodid tick Rhipicephalus appendiculatus transmits the apicomplexan protozoan parasite Theileria parva, which causes East coast fever (ECF), the most economically important cattle disease in eastern and southern Africa. Recent analysis of micro- and minisatellite markers showed an absence of geographical and host-associated genetic sub-structuring amongst field populations of R. appendiculatus in Kenya. To assess further the phylogenetic relationships between field and laboratory R. appendiculatus tick isolates, this study examined sequence variations at two mitochondrial genes, cytochrome c oxidase subunit I (COI) and 12S ribosomal RNA (rRNA), and the nuclear encoded ribosomal internal transcribed spacer 2 (ITS2) of the rRNA gene, respectively. RESULTS: The analysis of 332 COI sequences revealed 30 polymorphic sites, which defined 28 haplotypes that were separated into two distinct haplogroups (A and B). Inclusion of previously published haplotypes in our analysis revealed a high degree of phylogenetic complexity never reported before in haplogroup A. Neither haplogroup however, showed any clustering pattern related to either the geographical sampling location, the type of tick sampled (laboratory stocks vs field populations) or the mammalian host species. This finding was supported by the results obtained from the analysis of 12S rDNA sequences. Analysis of molecular variance (AMOVA) indicated that 90.8 % of the total genetic variation was explained by the two haplogroups, providing further support for their genetic divergence. These results were, however, not replicated by the nuclear transcribed ITS2 sequences likely because of recombination between the nuclear genomes maintaining a high level of genetic sequence conservation. CONCLUSIONS: COI and 12S rDNA are better markers than ITS2 for studying intraspecific diversity. Based on these genes, two major genetic groups of R. appendiculatus that have gone through a demographic expansion exist in Kenya. The two groups show no phylogeographic structure or correlation with the type of host species from which the ticks were collected, nor to the evolutionary and breeding history of the species. The two lineages may have a wide geographic distribution range in eastern and southern Africa. The findings of this study may have implications for the spread and control of R. appendiculatus, and indirectly, on the transmission dynamics of ECF.

Multi-locus genotyping reveals absence of genetic structure in field populations of the brown ear tick (Rhipicephalus appendiculatus) in Kenya.

Rhipicephalus appendiculatus is an important tick vector of several pathogens and parasitizes domestic and wild animals across eastern and southern Africa. However, its inherent genetic variation and population structure is poorly understood. To investigate whether mammalian host species, geographic separation and resulting reproductive isolation, or a combination of these, define the genetic structure of R. appendiculatus, we analyzed multi-locus genotype data from 392 individuals from 10 geographic locations in Kenya generated in an earlier study. These ticks were associated with three types of mammalian host situations; (1) cattle grazing systems, (2) cattle and wildlife co-grazing systems (3) wildlife grazing systems without livestock. We also analyzed data from 460 individuals from 10 populations maintained as closed laboratory stocks and 117 individuals from five other species in the genus Rhipicephalus. The pattern of genotypes observed indicated low levels of genetic differentiation between the ten field populations (FST=0.014±0.002) and a lack of genetic divergence corresponding to the degree of separation of the geographic sampling locations. There was also no clear association of particular tick genotypes with specific host species. This is consistent with tick dispersal over large geographic ranges and lack of host specificity. In contrast, the 10 laboratory populations (FST=0.248±0.015) and the five other species of Rhipicephalus (FST=0.368±0.032) were strongly differentiated into distinct genetic groups. Some laboratory bred populations diverged markedly from their field counterparts in spite of originally being sampled from the same geographic locations. Our results demonstrate a lack of defined population genetic differentiation in field populations of the generalist R. appendiculatus in Kenya, which may be a result of the frequent anthropogenic movement of livestock and mobility of its several wildlife hosts between different locations.